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Abstract 


An electromagnetic pulse (EMP) model is developed which allows a 


quick computation of the time development of the electric fields qene- 


rated by a high altitude nuclear burst. The model Is based on the 


Karza .-Latter high frequency approximation for high altitude EMP, which 


describes fields generated by Compton electrons interacting with the 


earth's magnetic field. The proper choice of a gamma time output func- 


tion, which can be integrated In closed form, and a small angle approxi- 


mation, made in the expressions for the Compton currents and air conduc- 


tivity, eliminate the time consuming numerical Integration% usually 


necessary in EP models to compute the Compton currents and air con- 


ductivity. This results in a considerable savings in computation time. 


The modes is presented in a manner which is simple to use but still 


allows the var 
'1on of the major theoretical parameters In the prob- 


lem. 


A simplified model of electron collision frequency as a function 


of electric field strength is given which enables the model to predict 


accurate results for nucleat weapon gamma yields up to at least 100 Kt. 


The results predicted by this EMP model compare to within 5.5% with 


results from the Air Force Weapons Laboratory CHEMP computer code. 


The computation time using the presented model on a CDC 6600 compu- 


ter is typically 5 sec or less for a 5 shake computation period in 


steps of .1 shake. 


The model presented should be useful for both classroom Instruc- 


tion and nuclear vulnerability/survivability studies and analysis prob- 


lems. 
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1. Theog 


Overview 


Since the solution to the EMP problem is actually the solution of 


a classical electromagnetic theory problem, the derivation of the moJel 


equations reduces to putting Miaxwell's equations into a convenient form. 


This is essentially accomplished by expressing Maxwell's equations in 


spherical coordinates and transforming to a retarded time frame. One 


must also develop expressions for the currents and conductivities of 


the system in the absorption region. The general equations dL.rlbing 


the high altitude model of Karzas and Latter have been derived in great 


detail by Chapman (Ref 3). 
Only the major points of the derivation will 


be given here. The system origin is assumed to be at the burst point 


with detonation at time t 
0 
0. This geometry is illustrated in Fig. 1. 


The key points to be remembered in this model are: 


1. Each gamma ray gives rise to one downward traveling 


Compton electron. 


2. The electrons are turned by the earth's magnetic field 


giving rise to a centrifugal acceleration. 


3. The relativistic electrons radiate energy in their 


forward direction. 


4. The gamma rays and the EMP radiation travel at the 


same speed. This leads to constructive interference 


of the radiation from each of the electrons. 


Particle Densities 


The gamma rays from 
ne nuclear weapon travel in a straight line 


to a point where they produce Compton electron,. 
At any given point r 


- 
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tie number of gammas which interact to produce Compton electrons Is 


pr. 
- 
dr] 
'TXr 
7 . 


E r= r2 )(r) 


where ' (r) is the mean free path of gamnia rays to produce Compton elec- 


trons, Y is the ga.na yield of the weapon in electron volts (eV), and E 


is the mean gamma energy in eV, 


Equation (1) may also be called the radial distribution function, 


or an attenuation function for interacting gaon'?a 
rays. The -term is 
E 


the total number of gar.,a rays available from the weapon. The °ar 
2 


term accounts for the divergence of the ga.ma rays as the radius r Is 


increased while the remaining terms account for the reduction in gaimmas 


due to th'.ir absorption in the atmosphere, based on the mean free path. 


It is assumed that the gala mean free path varies as the exponen- 


tial atmosphere. This gives the function3l relationship between A and 


r,- 


A(r) = A, exp [(HOB - r cos A).-, 
(2) 


where 


= gamma mean free path at standard pressure 


HOG 
height of burst in Km above the earth's surface 


r = radial distance from the burst Point to the point of interest 


A - angle between the position vector r and the vertical 


S - atmospheric scale height 


With this assumption, Eq (1) can be integrat.ed and becomes 


Y 
I 
S 
HODN 
r cos A,\ 
9(r -xp 
xp 
xp 
-S 
gr) = E 
Xr 
2 
exp 
A 0 cos A 
S. 
J 
' 


(3) 
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NJow if f(t) is the time distribution function of the weapon yield, 


the rate of Compton electrons, nc, produced at a given point r and time 


t is given by 


d nc u (r ) f t 
( ) 


Each Compton electron produces through inelastic scattering events 


several secondary electrons which forn the basis for the conductivity 


of the atmo" %ere. As in the Karzas-Latter approach, each Compton elec- 


tron is assumed to have a constant speed, V., throughout the ranne, R, 


of the electron which is a function of altitude. This allows the life- 


time to be expressed as R/VO . 
If each Compton electron produces secon- 


dary electrons at a constant rate, the rate of tecondary electron, ns, 


production is 


dns 
Ec/33 eV nc 


dt 
R/V-- nc 
(5) 


where Ec is the energy of the Compton electron and 33 eV is the average 


ionization energy per air molecule (Ref 2). 


Considering the differential current produced by the Compton elec- 


trons, it can be shown (Ref 3) that the Compton current and the number 


of Compton electrons are given by 


R.V0 
_ 
x (- N 
Jc,, 
= _ e g(r) 
v(z') F 
- 
+ 
dt 
(6) 


0 


R/Vo 
X 
(7) 
nc')=g(r) f 
f 
+ x (.r+ 
dT 


7 


.------ 
......... 


r.LP/Pli/,75- 13 


whecre 


r 


- 
the time since the creation of the Compton electron 


XW) - the radial distance the Compton electron has traveled 


e -the ma~gnitude of the electron charge 


The cuantity T is generally known as retarded time. 


It then follows from Eqs (4) and (7) that the number of secondary 


electrons is 


nFT 
-2Og(r) f 
f (T' 
!( 
dT~ 
dT 
(8) 


R- 
0 


where q is Ec/33 eV. 


-' 
Currents and Cor~ductivity 


In the Karzas-Latter theory, the s.peed of the Co-npton electrons is 


considered to be a constant, however, there is an acceleration due to 


the geo'-agnetic field. 
The general eiuat'on of notion for an electron 


in this case is 


d , mV=-e (E + VX B) meivcV 
(9) 
dt 


where 


=the electron rest mass 


V = the electron velocity 


E -the 
electric field 


B - 
the magnetic field 


vc the electron collision frequency 


Y = - 
(v0/c)2) 
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if the relativistic m~otion of th'p cle, tron is considcree, only the 


V A( 
8 term Is imnortant and in sol'crical coordinates, the expressons 


for the velocity co-nponents become (Ref 3) 


Vr 
V, (sin 2 9 Cos WT + cas2 0) 
(10) 


Ve 
V0 (COS e Sin e COS C.,- sin 6 cos o) 
(11) 


V0 - V0 (sin 0 sin ti) 
(12) 


where w I- the cyclotron frequency for an electron and is given by 


el0 
(13) 


~it 
athe mannitude of the geom1agnetic field. 


From Eq (10), X(-V) is found to be 


X (r) V0 sin 2 0 iB T 
+ vcos2 8e> 
(14) 


0) 
W 


The Compton currents may nowj 
be written as 


R/V0 


Jc(l) 
-eg(r) V0 I 
[f(T)(cos 2 0 + Sin e COS cWr)1 
d-r 
(15) 


r 
0 


RIN 


J' (T) - - e9(r) V0 f 
[f(T) sin 0 cos 0 (cos wV - 1)) 
dx' 


6 
0 


(16) 


R/V0 
Jc( ,) - 
- 
eg(r) V0 f 
[f(T) sin 0 Sin WTrI dT' 
(17) 


0 
0 


where 
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V° 
__sin 
____0 


T 
1 - 
T-- 
cos 
2 " 
+ 
v ° in2 
sin wT) 


In a sinilar nanner Eq (8) becomes 


q/v0 
- 
!ns( 
Vo 9(r) J Ff 
f(T-) dT-" 
d, 
(19) 


R 
.. 
o 
- 


where 


V" 
l- 
I --- cos2 0\ T"+V- 
sin2 0 
-- 
(20) 


C 
c 


: 
Equations (15), (16), (17), 
(18), and (20) may be simplified if 


the factor uT Is assumed to be small. 
In this case a Taylor series 


expansion for the sin and cos of aT, including only first and second 


order terms, is 


sin uT = (a 
(21) 


cos WT = I - 
-T 
(22) 


2 


The expressions for the Compton currents now become 


RV 
s 
2 
o 
dR/VO 


Jc(r) 
- 
ea(r) V0 J 
f(T) di' - 
sin 2 0 - 
r'2f(T) dT' 


0 
2 0 


(23) 


!2 R/Vo 


J 
) 
e(r) Vo sin e cos e 2--0 
/ 
V'2 f(T) dT" 
(24) 


R/V 0 


JC(T) - 
- eg(r) V0 sin 0 w j 
RV 
f(T) dT" 
(25) 


0 
0 


T -T- 
(I - 
0) ' 
(26) 
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where B = 
- 


In a like manner Eq (20) becomes 


T 
- 
(I - 
T" 
(27) 


An expression for the conductivity may also be found by using the 


equation of notion for the secondary electrons. These electrons are In 


the thermal regions iith energies ranging from about 10-15 eV to the 


ar.moient 
energy. 
It should be remembered for later use that the ambient 


energy of the secondary electrons is dependent on the electric field 


present. 
For consideration here, it is assumed that y -k I and also that 


the change of velocity with time is small compared to the other terms 


in Ea (9) so that d.%. may be neglected. Also, with low velocities, the 


dt 


V X 8 term is small compared to the remaining terms and may also be 


neglected. Then the velocity of the secondary electrons is 


V 
E 
(28) 


mVc 


Using Eq (28), the current due to the secondary electrons Is 


e2 
(2 


JS(T) =- e V ns(r) = 
E ns(r) 
(29) 
*c 


Co-lparing Eq (29) to JS = oE, an expression for the conductivity 


is 


e 
2 


o(r) - 
ns(r) 
(30) 
mv, 


Equations (19), (23), (24), (25), (26). (27), anl (30) provide the 


desired expressions for the Compton currents and the conductivity. 
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} 
Field Equations 


Maxwell's equations is rationalized UlS units are 


, 
at 


; 
" 
I aE 
-V x 
= UP + 
- 
(32) 


qv 
7.E = -- 
(33) 
£0 


V ' = 0 
(34) 


where q. is the total charge density and J is the total current density. 


In addition to these equations the continuity of charge requires that 


-qv 
+ v.; = 0 
(35) 


at 


Combining these equations to separate E and B and transforming 


them into spherical coordinates and into the retarded time frame (Ref 


3) the relations for E and B become 


" V2E+Ur 
1_.J 
. 
+ -- _Vqv 
r cc0 
£0 
(6 
(36) 


D -2 1 a 
+ 
0 
r (E) 
r)~ 
+ ~ 
~ 
- 
al' 
r 


2_ . 
a 
V 
j+ 
I 
.. 
rr (rB) 


(37) 


+_0 (uLJ 
0 - 
0 Joe)J =0 


In the Karzas-Latter model, only the time derivative portion of 


Eqs (36) and (37) are kept since the current variation with distance 


is slow compared to the variation in time for the high frequency com- 
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ponents. Also the fields an,! 
currents vary rapidly in time. 
This 


approxination is valid for about 100 shakes. Th s same hih freauency 


anoroximation is used here. In addition, the radial component of the 


field is dropped since it is weak compared to the transverse components 


and contributes only a very low frequency siqnal (Ref 2). 
The equa- 


tions for the transverse coinonents are 


9 
3 
21 
frEe 
+u 
Je,1 
0 
(8) 
0 


11- 
c(rB 
j 
. 0 
(39) 


-c r Dr 
c 
0 1 


2 
2 i 
21 . 
(rBe) +-.±- JO 
- 0 
(40) 


T1 
'Trr 
c 
J 


The currents in Eqs (38). (39), and (40) are total currents. The 


total currents are given by 


, 
.c 
+ a(T) 
(41) 
Oe, e, 
E 00 


Substitution of Eq (41) into Eqs (3,), (39), and (40) and integra- 


tion over time gives 


".12- 
(rEa + P 0 + VO O(T) 
Eq - 0 
(142) 


2i1 
C + PO O(T) E = 0 
(43) 


213 
c 
0 - U- (t) E - 0 
(4h) 


cr-E 
- 
0c 
Z' 


21 
(rB) 
+..J 
+ -% a(t) Ee -0 
(-5) 


cTr r 
+L 
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qutos(2)and (43) are In a forq wjhich can bc solved. 
The 


'1 
te"!' needed for solution of these enuat ions are the air conductivity 


-and 
the transverse components of the Compton currents. 
The air con- 


-~ductivity 
may be found by using Eqs (19) and (30). 
The transverse comn- 


nonents of the Compton currents may be found by using Eqs (24) 
and (25). 
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I'!- 
TYPICAL NARROW PULSE 
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IV. Resul ts 


To test the validity of the 1i,EMP m~odel, the ccwenuted values of 


this model were compared to available equivalent rodel values from the 


AFWL code for EMP calculations knowr as CIIE1P (Ref 14). 


The basic set o' conditions used for these calculations was: 


target location 
- nround zero 


height of burst 
= 100 Km 


geomagnetic field 
.3- 
nauss or 3 (lOP5 wb/rn 
2 


inclination asjvle 
= 0.6 degree 


Compten electron recoil energy - .75 MeV 


The oam.-a yield was vo, ?ed and the number of steps taken for the 


numserical integration wjas varied according to the ganxia yield, with 


mrore steps taken for the higher yields. 
Other pareameters varied for 


examination of Peal, field values were burst heirght, acmagnetic field, 


and pulse shape. A preionization level was also considered. 


The available data from the CHEMP(HI) code, which is CHEMP run with 


non-elf-consistent calculations, was computed using a pulse of the 


form of Eq (46) with a - 10- 
and 6 = 107 and the same geometry as 


given above. This pulse shane is almost identical to that of Eq O4R) 


with a - 107 and 6 = 3.7 (10)8. 
See Fig. 2 on Page 18. Usina this 


pulse shape, a range of gammva yields from . 01 Kt to 100 Kt was used to 


calculate the EMP field values. 
The peak field values are plotted in 


Fiq. 4. The values taken from the CHEMP(H) code are annotated by X 


marks. These peak field values show a maximum difference of 5.5% 


around the I Kt case, and a differzrnce of less than 2% for all other 


known cases. 
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